Clock shift in a strongly interacting two-dimensional Fermi gas 
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We derive universal relations for the radio-frequency (rf ) spectroscopy of a two-dimensional Fermi 
gas consisting of two spin states interacting through an S-wave scattering length. The rf tran- 
sition rate has a high-frequency tail that is proportional to the contact and displays logarithmic 
scaling violations, decreasing asymptotically like 1 / {uj'^ In^ to) . Its coefficient is proportional to 
ln^(a2£)/a2D), where 020 and a2j5 are the 2-dimensional scattering lengths associated with initial- 
state and final-state interactions. The clock shift is proportional to the contact and to In(a2£)/a2i5). 
If I ln(a2£)/fl2D)| 2> 1, the clock shift arises as a cancellation between much larger contributions 
proportional to In^ {a2jy /a2D) from bound-bound and bound-free rf transitions. 
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Spectroscopy of the internal levels of atoms is a cen- 
tral subject of atomic physics. The associated transition 
frequencies of a single atom are known with the precision 
of atomic clocks. In fact, atomic clocks use a hyperfine 
transition in cesium to keep our standard of time. In a 
gas, the atoms experience energy shifts due to interac- 
tions with the surrounding atoms, resulting in transition 
frequencies that differ from those of free atoms. Such 
density-dependent interaction shifts are a major source 
of systematic errors in atomic clocks, and are thus called 
clock shifts. Precision spectroscopy aims to avoid these 
shifts as far as possible, e.g. by working at extremely 
low densities or by using a purely spin-polarized gas of 
fermions in the ultracold limit where the remaining S- 
wave collisions are forbidden by the Pauli principle [1]. 
Alternatively, the atoms may be stored in an optical lat- 
tice, where both the center-of-mass motion and the in- 
teractions are completely quenched [5]. 

In investigations of many-body physics using ultracold 
gases, the clock shift is, however, a signal of interest 
rather than a nuisance. It can be especially interest- 
ing when the atoms are strongly interacting, i.e. their 
scattering length is large compared to both the range of 
interactions and the average interparticle distance. The 
interaction energy Jj i3j and the pairing gap [1] have been 
measured for strongly-interacting fermionic atoms near a 
Feshbach resonance. The generic setup for these experi- 
ments involves a gas containing atoms in just two hyper- 
fine states |1) and |2). A radio-frequency (rf) pulse, tuned 
to the hyperfine splitting between |2) and a third, unoc- 
cupied state |3) will then transfer atoms from |2) into |3). 
If, for example, states |1) and |2) form a molecular bound 
state, the rf pulse has to be detuned compared to the bare 
hyperfine splitting by a frequency uj associated with the 
binding energy of the molecule. At the two-body level, 
effects like these can be analyzed in a straightforward 
manner, even in the presence of strong final-state inter- 
actions between |1) and |3) [5]. In a genuine many-body 
situation, however, analytic methods for calculating the 



rf transition rate r(w) are not available. Exact relations 
for r(aj) that hold independent of density and interaction 
strength are therefore of considerable interest. 

For strongly interacting fermions with two spin states 
in three dimensions (3D), a number of exact rela- 
tions that connect thermodynamic variables to large- 
momentum and high-frequency tails of correlation func- 
tions have been derived by Tan [6l [7]. These relations 
all involve the contact C, which basically measures the 
probability for pairs of fermions in states |1) and |2) to 
be very close together. More precisely, the contact may 
be defined by the universal power-law tail 

n^(fc)^C/fc^ (7 = 1,2 (1) 

of the momentum distribution of either spin state at large 
momentum k. The Tan relations are universal in the 
sense that they apply to any state of the system, e.g. 
few-body or many-body, homogeneous or in a trapping 
potential, Fermi-liquid or supcrfluid state, provided only 
that the length scales associated with the temperature 
and number densities are large compared to the range of 
interactions. The origin of this universality is that the 
Tan relations are a consequence of operator identities, 
some of which follow from the operator product expan- 
sion of quantum field theory [S]. There are also uni- 
versal relations governing rf spectroscopy. Specifically, 
the average clock shift, i.e. the first moment of r(a;), is 
proportional to the contact C [9] or, equivalently, to a 
derivative of the total energy with respect to the scat- 
tering length [TD]. Moreover, the rf transition rate r(a;) 
has a high-frequency tail that is proportional to C and 
decreases like 1/w", where the exponent n is either | or 
I depending on the strength of the final-state interac- 
tions between |1) and |3) [TTHI3] . Some of the universal 
relations have been verified experimentally, e.g. by com- 
paring the values of C obtained from the tail of the mo- 
mentum distribution in Eq. ([T]) and from the rf transition 
rate at large w [Ti] . 

In our present work, we derive universal relations for 
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rf spectroscopy in fwo-dimensional Fermi gases. These 
systems have been investigated in several recent experi- 
ments [I5HT7] . Beyond the motivation provided by these 
experiments, universal relations that constrain the rf 
transition rate in two dimensions (2D) are of interest 
also for other reasons. First, rf spectroscopy measure- 
ments can be performed within a rather short time scale. 
It thus allows access to not only equilibrium but also ex- 
cited states of the many-body system, e.g. the polaron on 
the repulsive branch of a Feshbach resonance [IB] or Bose 
gases in the regime of strong repulsion. A second reason 
is that a 2D gas with zero-range interactions provides an 
example of a non-relativistic many-body system that is 
scale invariant at the classical but not at the quantum 
level. This breaking of the scale invariance by quantum 
effects is an elementary example of an anomaly [19]. It 
implies that the coupling constant g (which will be de- 
fined explicitly below) is in fact a running coupling con- 
stant .g(/i) = — 27r/ ln(a2£)/^) that changes logarithmically 
with the momentum scale /x, reminiscent of the coupling 
constant of Quantum Chromodynamics. The associated 
intrinsic length a2D sets the scale for the two-body scat- 
tering amplitude which, at low energies, has the charac- 
teristic form EOl El] 



/(?) 
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(2) 



For any positive value of 02,0, this amplitude has a pole 
at g = i/^2D, which signifies the presence of a diatomic 
molecule with binding energy Efi = h^/ma^j^ that we 
will call the (12)-dimer. 

In experiments with ultracold atoms, the 2D gas arises 
from a 3D system by adding a strong transverse confin- 
ing potential that restricts the atoms to the ground state 
in the third dimension. For harmonic trapping with an- 
gular frequency , the relation between ^nd the ex- 
perimentally tunable 3D scattering length a is given by 
£/a — fi{i'^ / 0,20)1 where = h/muJz and the function 
fi{x) is given in section V.A of Ref. f21i . 

A system consisting of fcrmions with the two spin 
states |1) and |2) can be described by a quantum field 
theory with the interaction Hamiltonian 



Hint = ig/m) J d^R i^lt^y2MR 



(3) 



n„{k) behaves asymptotically like C /k^ with no logarith- 
mic corrections, as in 3D and also in ID [22]. The contact 
C — J (P R C (R) is a.n integral of the contact density 



C{R)^g'{44ij2MR)), 



(4) 



which is the expectation value of a local operator that an- 
nihilates and recreates a pair of fermions. The structure 
of these results is identical to the 3D case and, indeed, 
many of Tan's universal relations may be generalized to 
2D [23 . Universal relations for current correlators, struc- 
ture factors, and hydrodynamic correlators in 2D have 
also been derived using the OPE [Jf . 

We now consider the rf spectroscopy of the 2D system 
consisting of fermions in states |1) and |2). Beyond the 
interaction described by 020, we also include final-state 
interactions between |1) and |3) with 2D scattering length 
'^2D ■ will refer to the diatomic molecule consisting of 
|1) and |3), which has binding energy E'^ = fi? / ma^jj'^ , 
as the (13)-dimer. The transfer rate r(a;) between |2) 
and 1 3) can be expressed in terms of the imaginary part 
of a time-ordered correlation function: 

where is the Rabi frequency of the rf transition and 
0^i{uj, R) is an operator that depends on a complex vari- 
able uj: 



xT4i;3{R+^r,t)ijl^2{R-lr,0). (6) 

(The time-ordering symbol T acts on all operators to the 
right.) The bilocal operator can be expanded in terms of 
local operators by using the OPE: 



0,.t{LJ,R) = j:^Cn{io)OJR). 



(7) 



Here g is the bare coupling constant, which can be iden- 
tified with the running coupling constant at the momen- 
tum scale of the ultraviolet cutoff. Within this frame- 
work, universal relations for both equilibrium properties 
and rf spectra can be derived using the operator prod- 
uct expansion (OPE) [S]. In particular, the tail of the 
momentum distribution in Eq. ([T]) follows from the OPE 

hi tpl{R+ lr)i;„{R- ^r). The leading contribution {0,i{Lj,R)) > -n2{R) 



Local operators can be assigned scaling dimensions. The 
quantum fields ipa- have dimension 1. The gradient V and 
the time derivative d/dt increase the dimension by 1 and 
2, respectively. At large frequency w, higher dimension 
operators in the OPE have Wilson coefficients C„ (lu) that 
decrease with higher powers of ui. The Wilson coefficients 
Cn{u}) in Eq. ^ can be calculated by matching Green 
functions of the bilocal rf operator on the left side of 
the OPE with Green functions of the local operators on 
the right side. For operators of dimensions 2 and 4, it is 
sufhcient to match the Green functions in the 1-atom and 
2-atom sectors. Our final result for the expectation value 
of C'rf(w,i2) at large complex w, including all operators 
with dimensions up to 4, is 



i V2n2(i?) +4mV • J2{R) 



that is non-analytic as r — )■ comes from the operator 
V'i'02V'2V'i(-R)j whose coefficient exhibits a logarithmic 
singularity in 2D of the form |rplog|r|. As a result. 



+ - 



iln(a2£,/a2_D) \n{a2DV—rmj , 



C{R), (8) 
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where n2 and J2 are the number density and current 
density for |2). If the system is in a steady state, V • J2 
vanishes because of current conservation. The integral 
of V^n2 vanishes by the divergence theorem. The terms 
shown exphcitly in Eq. (|8| are not affected by final-state 
interactions between |2) and |3), which first enter into 
the coefficients of dimension-6 operators. 

We can derive various universal relations from the OPE 
in Eq. ([s]). The tail of the rf transition rate at large 
positive frequencies w can be obtained by extrapolating 
Eq. ([8| to the real axis, taking the imaginary part, and 
then inserting it into Eq. ([s]): 



4mcj2[ln2(w/£;^) -|-7r2 



(9) 



Higher dimension operators in the OPE give contribu- 
tions to the high-frequency tail that decrease faster than 
Lu^^ . These contributions are suppressed by a posi- 
tive power of an energy scale of the many-body sys- 
tem divided by u). The relevant energy scales of the 
many-body system include Ed, the thermal energy ksT, 
and the two Fermi energies ^n'^ni/m and ATT^n2/m. 
The factor \n^{E'jEd) in Eq. H can be interpreted as 
VoT^^ I (\ I OLr^Y)) , where ^(l/ajjjy is the running coupling 
constant at the momentum scale l/a^jj- If 0.20 ~ ^2D, 
the high-frequency tail in Eq. ^ vanishes identically. 
In this case, r(w) just consists of an unshifted peak 
7rri2iV2(5(w), because the rf pulse merely rotates |2) and 
1 3) in spin space [TU]. The nontrivial dependence on uj 
m Eq. (|9| is a consequence of the anomalous scale in- 
variancc. If |ln(a;/i?^)| ^ tt, the high-frequency tail in 
Eq. ^ has the simple scaling behavior r(w) ~ l/w^. 
However at much larger or much smaller frequencies, 
there are logarithmic scaling violations and the asymp- 
totic behavior is T{w) ~ {w"^ It? w)~^ . Simple scaling 
behavior also arises in the limit of negligible final-state 
interactions in which ln(i?^/aj) 00. In this limit, 
the high-frequency tail in Eq. ^ reduces to fl'^C/Amuj'^. 
This relation has been used in recent experiments with 
""^K atoms to extract the contact density in the normal 
state of the 2D Fermi gas [3S] . The contact density and 
its dependence on 0213 is close to that obtained from vari- 
ational Monte Carlo calculations at T = [26] . 

Sum rules can be derived by expressing weighted inte- 
grals of r(w) as contour integrals in the complex-w plane: 



+ C30 

du! f {uj)r {oj) 



dPRj>dujf{uj){0,i{uj,R)), (10) 



where the u contour runs from +00 to —00 just below 
the real axis and then from —00 to +00 just above the 
axis. If the contour is deformed into a circle of infinite 
radius, the OPE in Eq. (|8| will be exact everywhere on 
the contour. Choosing f{uj) = 1 and /(w) = w, we obtain 




- 8 2 4 6 

FIG. 1. Rf transition rate for the (12)-dimer as a function 
of u)/Ed for the case E'^ = 9.4i5d. There is a delta function 
a,t u) = —{E'd — Ed) from the bound-bound transition. The 
threshold for bound- free transitions is ai lj — Ed- The grey 
dots in the bound-free region are data from Zwierlein and 
collaborators [27], normalized by a least-squares fit to the 
theory curve. 
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(lib) 



These sum rules are exact: they receive no contributions 
from higher dimension operators in the OPE. The first 
sum rule guarantees that interactions do not change the 
integrated rf signal. The second sum rule determines 
the clock shift (w), which is defined as the ratio of the 



integrals in Eqs. (Ill: 



HE'jEd)C 



(12) 



Note that both the clock shift and the tail of the rf tran- 
sition rate in Eq. ([9| are unaffected by final-state inter- 
actions between |2) and |3). 

The universal relations in Eqs. (|9| and (11) can be il- 
lustrated by considering rf spectroscopy on the two-body 
system consisting of a single (12)-dimer with binding en- 
ergy Ed- The rf transition rate is 



r(c.) = 



Trn^E'^Edln'iE'jEd) 



{E'd-Ea? 
nQ^Edln\EyEd) 
u^[ln^{{u-Ed)/E'J + 7:^ 



Siu; -Ed + E'd) 

eiuj-Ed)- (13) 



By evaluating the integrals in Eqs. (11), we verify that 



the sum rules are satisfied exactly with iV2 = 1 and C = 
4:TrmEd, the contact for the (12)-dimer obtained from the 
adiabatic universal relation dEd/da2D = — C/(27rma2D) 
[23!. The universal relation for the tail in Eq. ^ holds 
for |a;| ^ Ed- The clock shift for the (12)-dimer has a 
negative contribution proportional to \n^{E'^/Ed) from 
bound-bound transitions. If \ \n{E'j^/Ed)\ ^ 1, there is 
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a large canceling contribution from bound-free transi- 



tions, such that the total clock shift agrees with Eq. ( 12 ) 



Near the threshold, the bound-free transition rate has 
the behavior l/ln^((a; - Ej)/E'j). This follows from the 
inverse-logarithmic behavior of the two-body scattering 
amplitude for |1) and |3) analogous to Eq. Recent 
experiments using ^Li atoms in 2D II 7j are in good agree- 
ment with the bound-free transition rate in Eq. ( |13[ ) for 
a single (12)-dimer in the presence of final-state interac- 
tions. The data in Fig. [l] [23 were measured near the 
3D Feshbach resonance for |1) and |2) at a magnetic field 
of 690.7 G and a lattice depth to recoil energy ratio of 
18.6(7), and they extend out to higher u) than in Ref. [TT] . 
The ratio E'^/Ed is around 9.4 and Ed is about a factor 
6 larger than the Fermi energy [37]. This places the sys- 
tem in the BEC limit of the many-body problem where 



the two-body result in Eq. ( 13 ) applies 



We now consider the many-body system in the limit in 
which E'^ is much greater than Ed and the other relevant 
energy scales. The bound-bound and free-bound region 
is then well separated from the free-free and free-bound 
region of the rf spectrum. There could be additional 
structure in the region of w between —{E'^ — Ed) and 
associated with clusters of three or more atoms. However, 
if this region is featureless, it is useful to define a clock 
shift for the rf transition rate restricted to a range of 
frequencies — wo<w<i^o that excludes the bound-bound 
peak. We choose ujq and E'^ — ujq to be much greater than 
the energy scales of the system. This clock shift is given 



by a ratio of integrals analogous to those in Eq. (Ill 



They can be expressed as contour integrals analogous to 
those in Eq.(lO), except that the contour runs from +ujq 



to —ujQ just below the real axis and then back to +a;o just 
above the axis. If the contour is deformed into a circle of 
radius ojo, we can apply the OPE in Eq. ([8| and obtain 
the sum rules 



+1^0 

du! r(w) 



+1^0 

du ujT{uj) 



In^ 

a2D 



In^ 
a2D 



ImEi(-Lo 
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TTmE',, 



TT 



"'2D 
a2D 



arctan 



in 



(14a) 



2m 
(14b) 



where Ei(a;) is the exponential integral function and 
Lq = ln(cJo/£'d)- C)ne advantage of finite- frequency sum 



rules such as those in Eqs. ( 14 1 is that they can be evalu- 



ated using measured rf spectra without extrapolations to 
larger frequencies. As indicated by ..." in Eqs. (14), 



these sum rules are not exact. They receive contributions 
from higher dimensional operators in the OPE that are 
suppressed by powers of an energy scale of the system 
divided by loq or E'^^. For example, the higher dimen- 
sion contribution in Eq. ( 14a I is suppressed by a fac- 
tor of C jmE'^Ni. The corrections that are suppressed 



by inverse powers of wq can be minimized by choosing 
as large as possible while still avoiding contributions 
from the bound-bound peak. Neglecting the contribu- 
tions from higher dimension operators, the clock shift 



obtained from Eqs. (14) reduces to 



In 



2 ^ 
Ed 



arctan 



\U\ 



-In^^ 



C 



Ed I 47rmiV2 



(15) 



If \ \x\{E'J Ed)\ ^ 1, this clock shift is much larger than 



that in Eq. (12). The clock shift in Eq. (15), which ex- 



cludes the contribution from bound-bound transitions, is 
the one that is more relevant to the limit of weak final- 
state interactions. The factor hiiE'^jEd) diverges in the 
limit al^Y) 0, indicating that the clock shift in this case 
is quadratically logarithmically sensitive to the range of 
interactions. 

In summary, we have presented universal relations for 
the rf spectroscopy of strongly interacting fermions in 
a 2D Fermi gas. They are of direct relevance to cur- 
rent experiments with ultracold atoms, providing rigor- 
ous connections between rf spectrocopy, thermodynam- 
ics, and various correlation functions. They reveal as- 
pects of these many-body systems that are directly re- 
lated to the breaking of scale invariance in 2D and will 
hopefully provide inspiration for further investigations of 
these unusual systems. 
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